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Asymptotic  Theory  of  Second  Order  Differential 
'Equations  with  Two  Simple  Turning  Points 

Nicholas  D.  Kazarinoff 

Communicated  by  L.  Cesari 

i  Introduction1.  In  1931  hanger  initiatcrl  and  gave  the  first  of  numerous 
contributions  to  what  has  become  a  successful  theory  for  asymptouc  expanse™ 
of  the  solutions  of  a  differential  equation  with  a  turning  point.  This  theory  ), 
been  extended  and  appliS  to  a  great  many  questions  by  him  and  by  othc^ 
An  extensive  list  of  references  may  be  found  in  L .  Cesari  s 
behavior  and  stability  of  soSjions  of  differential  equations ,  Chap  e  ,  p  g 

Verlag,  Berlin,  1958.  \j  ht  ;■  r  ■  -  ”***  ' 

a  generaiasymptolrcrth^ryW  ordinary  second  order 

JZ  SfaSTSSorflth  ^Simple  turning  fan<s  and  conta,n,ng  a 

mimericallv  large  parameter.  IrrpartiCTfe^a^^  *  ^symptom 

the  form 

Wc  consider  this  equation  for .  in  a  dosed, 

region  9,  of  the  complex  plane.  We  assume  that  f«\X\>N  and  for  sQA 
coefficient  P{s,  X)  is  of  the  form  ■  \ 

•2  */(*>*”'. 


(1.2) 


P(s,  X) 


*hm  ecu  h  *  analytic,  and,  most  importantly,  that  *<»  **  <“ 

simple  zeros,  a  and  0,  in  the  interior  of  2 

In  a  region  which  includes  a  turning  point,  i.e..  a  zero  or  sini  ulantj  of  p0{s 
tire  solutions  of  the  differential  equation  (t.t)  depend  upon  ?.  m  so  tntmatc  a 

-  ,  Slates  Air  Force  through  the  Air 

P^s^iSirc^^  Com,n 

**  s""ic  symhota  ior  '”"uvc 

constants. 
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which  includes  just  one  simple  turning  point  are  known  and  are  based  upon 
Airy  functions  [2,  5j.  Thus,  it  would  appear  that  the  behavior  of  solutions  of 
(1.1)  ov(;r  the  entire  region  which  contains  two  simple  turning  points,  is 
obtainable  by  the  familiar  procedure  of  evaluating  the  coefficients  in  the  depend¬ 
ence  relations  which  connect  solutions  whose  behavior  is  known  about  one  of 
the  turning  points  with  those  solutions  whose  behavior  is  known  about  the  other 
turning  point.  Unfortunately,  this  evaluation  is  not  possible  in  general;  two 
linearly  independent  solutions  may  have  the  same  dominant  asymptotic  form, 
which  makes  the  inference  of  the  identity  of  two  solutions  from  the  identity  of 
their  asymptotic  forms  invalid.  Thus,  a  new  theory  is  necessary  if  we  desire  to 
have  uniform  asymptotic  expansions  of  solutions  of  (1.1)  over  all  Langer  [4] 
has  derived  the  leading  terms  of  such  asymptotic  expansions  in  the  special  case 
where  s  is  a  real  variable  on  a  bounded  interval,  p0 (s)  is  real-valued,  and  pl  (s)  =s0. 
We  derive  the  asymptotic  expansions  to  n  -4- 1  terms,  where  n  is  an  arbitrary 
non-negative  integer,  of  the  solutions  of  (1.1)  under  the  general  hypotheses  set 
forth  in  the  first,  paragraph  above  together  with  some  others  of  a  more  technical 
nature  to  be  set  forth  later. 

Interest  in  the  problem  discussed  here  stems  mainly  from  possible  applications 
for  the  theory  derived.  In  certain  regions,  the  differential  equations  for  the 
angular  and  radial  spheroidal  functions  are  of  the  type  (1.1).  This  is  also  true 
of  the  Whittaker  equation  for  certain  configurations  of  its  parameters.  The 
spheroidal  functions  are  important  in  problems  concerning  scattering  by  a  prolate 
spheroid.  The  Whittaker  functions,  disguised  as  Coulomb  wave  functions,  occur 
in  quantum  mechanics.  Equations  of  type  (1.1)  are  also  of  interest  in  other 
problems  of  wave  motion  and  diffraction. 

We  have  divided  the  discussion  below  into  three  parts.  In  Part  1  we  transform 
the  differential  equation  (1.1)  into  one  more  suitable  for  analysis.  We  call  this 
canonical  form  the  given  equation.  We  then  give  an  algorithm  for  the  construction 
of  a  related  equation  whose  coefficients  resemble  the  coefficients  of  the  given 
equation  to  an  arbitrarily  prescribed  degree.  In  Part  II  we  study  the  solutions 
of  the  related  equation.  These  involve  Weber  functions  of  large  complex  order 
and  argument.  They  have  been  studied  by  Erdklyi,  Kennedy  &  McGregor  [3]. 
We  make  considerable  use  of  their  results  and  have  shown  that  their  asymptotic 
representations  bold  uniformly  in  an;  v  over  a  finite  range  of  arg  v,  where  v  is 
the  order  of  the  Weber  functions  h  volved.  We  also  give  an  algorithm  for  re¬ 
cursively  determining  the  terms  in  the  asymptotic  expansions  of  these  Weber 
functions.  In  Part  III  wc  prove  that  solutions  of  the  related  equation  are  asympto¬ 
tic  expansions  to  n  -f  1  terms  of  solutions  of  the  given  equation,  n  being  an’- 
non-negative  integer.  The  method  of  proof  is,  as  usual,  to  transform  the  given 
equation  into  an  integral  equation  of  Volterra  type,  whose  kernel  involves  solu¬ 
tions  of  the  related  equation,  and  to  solve  this  integral  equation  by  the  Picard 
method  of  successive  iteration.  Wc  also  give  approximations  for  derivatives  of 
solutions  of  the  given  equation.  Our  main  results  are  stated  as  Theorems  1  and  2 
of  §§  12  and  1 ).  An  especially  interesting  feature  is  the  occurence  of  a  denumerable 
number  of  characteristic  solutions  of  the  given  equation.  These  are  bounded  and 
oscillatory  on  certain  curves  joining  the  turning  points  and  exponentially  small 
on  their  extensions  beyond  tire  turning  points. 
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Part  I:  Construction  of  the  Related  Equation 
2.  The  given  equation.  The  analysis  to  follow  involves  certain  functions 
and  a  mapping  connected  with  p0  (s) .  To  simplify  their  form  we  adopt  a  normalized 
form  of  the  differential  equation  (t.l)  as  the  basis  of  the  sequel.  This  normaliza¬ 
tion  was  used  by  Langes  in  [4,  §§  2,  )].  Its  adoption  here  will  Ik  seen  to  entail 
certain  assumptions  on  />0(s). 

We  consider  the  mapping  from  the  given  region  onto  a  region  of  the 
z-plane  defined  implicitly  by  the  equation 


2.1) 


/(**-  Wdt-  fpl{t)dt~0. 
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An  immediate  question  is  whether  or  not  (2.1)  defines  a  mapping  at  all,  and  if 
so,  is  it  a  schlicht  mapping?  We  can  shew  that  (2.1)  defines  a  mapping  and  that 
at  each  point  of  Qi  it  defines  a  locally  one-to-one  mapping.  The  function  which 
is  the  left  member  of  (2.1)  is  analytic  in  z  except  at  -J  1  and  in  s  except  at  a  and  ft. 
Its  partial  derivatives  do  not  varysh  except  on  the  lines  z  —  ±1  and  s~a  or  p. 
Thus  by  the  implicit  function  theorem  for  analytic  functions,  there  exist  in 
neighborhoods  of  all  points  of  &i%  except  possibly  a  and  ft,  analytic  solutions 
z(s)  of  (2.1)  with  inverses  s  (z).  Further,  d sjd z  (s*  —  1 ) 4  1  (s) ;  and  hence, 

s'(*)=F  0  except  perhaps  at  ±1. 

We  now  examine  what  happens  at  the  exceptional  points  a  and  ft.  A  com¬ 
putation  shows  that  we  may  write 


(2.2) 


/V 


i)^dt  =  (z  —  1)*$i(z)  and  //>*(/)  dt  =  (s  ~  a)3 1]  (s) , 


where  $>,  and  I]  arc  analytic  in  neighborhoods  of  -f- 1  and  a,  respectively,  and 
neither  <£,(1)  nor  Pi  {%)  is  zero  Thus,  we  may  write 


(2-3) 


F(z,  s )  =  (z  - 


-A 

1  \*i 


(*-*)  *?(*)  =  0, 


where  F  is  analytic  in  a  neighborhood  of  (1,#)  and  neither  dFjdz  nor  dFjds 
vanishes  at  (I ,  os).  Therefore,  by  the  implicit  function  theorem,  there  exists  a 
solution  z[s)  of  (2-3),  and  hence  of  (2.1),  with  inverse  s(z),  which  is  analytic  in 
a  neighborhood  of  a  and  such  that  z(tx)  —1  and  s'(l)=|=0. 

In  order  to  apply  this  same  technique  to  (2.1)  at  s—ft,  we  must  assume  a 

-10 

normalization  of  the  coefficient  p0.  The  integrals  /  {P—i^dt  and  /  f\{t)  dt 


are  independent  of  path  provided  the  paths  do  not  encircle  the  singularities  of 
the  integrands.  However,  it  may  be  that  the  latter  integral  vanishes.  We  assume 
the  contrary  and  choose  the  parameter  so  that* 

(i)  / p\{t)  dt  —  /  {P  -  i)idt  =  - 

a  1 

W'e  can  now  apply  considerations  analogous  to  those  above  to  reach  the  con¬ 
clusions  that  (2.1)  has  an  analytic  solution  z(s)  at  (--!,/?)  with  inverse  s( a 

s  We  do  not  consider  the  limiting  behavior  of  the  solutions  of  (l.i)  as  the  turning 
points  a  and  ft  approach  each  other.  F.  W.  J  Olver  has  informed  the  author  by  a 
private  communication  that  he  is  examining  this  limiting  behavior. 
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whose  derivative  s'  is  nonzero  at  ~  1.  Jf  arg  p0  is  constant  on  the  line  segment 
joining  a  ami  ft,  the  hyjHithesis  (i)  is  clearly  fulfilled.  This  is  the  case  treated 
by  Larger  [/];  ami  since,  in  this  instance,  the  mapping  defined  by  (2.1)  is  one 
to  one  on  the  line  segment  joining  a  to  fi,  it  is  schlicht  in  an  open  neighborhood 
of  this  closed  segment. 

We  have  shown  that  the  implicit  relation  (2.1)  defines  a  mapping  from  9>t 
onto  a  region  9S  and  that  i‘  defines  a  locally  one  to  one  mapping  of  some  neighbor 
hood  of  any  point  of  9t  onto  a  neighborhood  of  9..  We  require  more  than 
this  and  assume  that 

(ii)  The  relation  (2.1)  defines  a  schlicht  mapping 


;(s): 

with 

Under  this  assumption,  the  change  of  variables 


(2.4a) 


s  «  n{z) , 


(2.4  b) 


>*(*)  “  {-Sf }*«(*)• 


where  s  satisfies  (2.1),  transforms  the  differential  equation  (1.1)  into 


A*0M) «  =0. 


where  Q  is  analytic  for  j  A  J  >  N  and  for  z£9..  We  write 


In  particular, 


9o(«)  =(**-<)• 


We  henceforward  refer  to  the  differential  equation  (2.5)  as  "the  given  equation”.. 


3.  The  first  approximating  equation.  The  analysis  of  the  given  equation 
is  based  upon  the  construction  of  an  equation  which  resembles  it  up  to  terms  of 
the  form  A-n_!0(1)  in  the  coefficient  of  «,  where  n  is  any  non-negative  integer 
and  0(1)  denotes  a  function  of ;  and  A  which  is  bounded  for  2<~ 9>t  and  for  |  A|>iV. 
The  algorithm  lor  the  construction  of  this  related  equation  is  similar  to  the  algo¬ 
rithm  of  McKelvey  given  in  [6].  The  matter  at  issue  in  McKelvey’s  paper 
is  the  approximation  of  solutions  of  an  equation  with  a  single  turning  point  of 
order  two.  Such  a  turning  point  may  be  thought  of  as  the  confluence  of  two 
simple  turning  points  so  that  similarities  in  the  analyses  of  the  two  cases  are 
not  unexpected. 

The  point,  of  beginning  in  both  instances  is  Weber’s  equation,  which  may¬ 
be  written  in  the  form 


V  =  0. 


£2} V IV*  ( m  .  „ 


mii'iwxu  ’J»+*tYKiT:  ffax+nri *efiw.  - 1 


■u.  **.,  U  *<t,t«UV»jrtS*f'lri"»  w»tw  ■»<*•«! 


1- SfV^Tt  *■  **«».<- 


rf*r 

r/;1 
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We  make  the  transformation 

/“(2A)*2,  !  "(/}  «  »(*) , 

and  Weber’s  equation  becomes 

•A(A*2~r)r  ~  0. 

Let  a  non-negative  integer  n  be  chosen  once  for  all.  We  now  choose  v  so  that 

(3-1)  v  —  A  -*  c  —  2]  Cj \ 

0 

where  the  (arbitrary)  constants  c  and  r;  are  yet  to  be  determined.  This  brings 
Weber’s  equation  into  the  form 

(3-2)  0~[a4(s2-1)-FAc-4-£ c;-;.-']t’-0  or  f"-/?foIA)r-0, 

which  is  our  first  approximating  equation. 

4.  The  second  approssntaring  equation.  Wc  next  construct  an  equation 
resembling  the  given  equation  in  both  the  A2  and  A  terms  of  the  coefficient  of  u. 
Formally,  this  step  is  almost  identical  with  the  corresponding  one  in  [£?].  The 
idea  is  to  make  a  change  of  dependent  variable 


(4.1) 


Z  ~fi0v  v'lX 


in  the  first  approximating  equation  (3.2)  and  to  determine  the  functions  ft0 
and  /q  in  such  a  way  as  to  make  the  new  differential  equation  in  Z  more  closely 
resemble  the  given  equation.  By  differentiating  Z  twice,  employing  (3.2)  at  each 
step  to  replace  v"  by  Rv,  and  by  constructing  the  eliminant  for  v  and  v'  among 
the  three  relations  connecting  Z,  Z' ,  and  Z"  with  v  and  v' ,  we  find  that  the  dif¬ 
ferentia!  equation  satisfied  by  functions  Z  of  the  form  (41)  is 


(4.2) 

wherein 

(4.3  a) 

(4.3  b) 
and 
(4-3  c) 


d*Z 

dz* 


H  dZ 

dz 


1.Z  —  0 
D0  dz  D„ 


E>0(*,A)» 

H(z,a)  - 

/M)  = 


Fo  /<o  +  lh  #/A  j 

\fhl%  Fa+/q'M 

/q/A  2/4+ /q/A  j’ 

,«o  +  /q  R  +  2/q  R  A'1  +  /q  RjX  /< o  +  /q  /?/A 
2/i'o  +  ft"  A'1  +  /q  Rj/.  /q  +  /4M 


In  so  far  as  A  is  concerned,  the  functions  D0  and  H  arc  bounded  for  |A|>A' 
with  D0(c,  o.  -  •  n\q0. 

We  shall  1  rinine  /q  and  //,  by  the  condition  that  the  terms  in  the  coefficient 
y/D0  in  (4.2!  '  ich  become  infinite,  with  A  have  sum  A s?f.(*)  +  A?»(;:).  We  recall 

that  the  choic  4  the  constant  c  in  ( }.2)  is  in  our  hands.  Now,  a  simple  division 
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shows  that 


/  _  P  t  n-i  /'o  +  2/q  M-l+  /q  R'i>.  !>o  +  /q  A’M 
l>o  l"  *  2f<'0 +„"!?.  Ha+HiH  ' 


Recalling  the  definition  of  R  from  (J.2),  we  observe  that  the  coefficient  of  in 
/  Da  is  q0.  The  coefficient  of  /.  may  be  computed  and  is 


c  -f  I>o’  (z>  so) 


2/'i?o+/V/o  » i?o 


D^l(r,  oo)  [{c/4  +  2/q,/q  qQ  +  /'o/'i  ?o)  ~  (<7'i  7o  +  2/iJ/t,  7-0]  • 


Therefore,  a  sufficient  condition  that  this  coefficient  be  equal  to  </,  is  that  /q> 
and  /q  satisfy  the  following  system  of  differential  equations: 


2/q  <h  +  lh  <h  r=  /'o(?r 

2/4  =/'j{7i“c)- 


with  boundary  condition  that 


/<o—  =  °°)  =1- 

A  possible  solution  of  this  system  is 


(4.5a) 


//„  —  COSh  {}  ,  fly 


sinh  & 


where 


(4.5  b) 


ip(z)  =  (.:3  —  1)},  arg<jn  =  0  for  2>1, 


(4.5  c) 


{Hz)  ^  f$'ls)f£ds. 
'  J  2  <p(s) 


We  deduce  from  these  formulas  that  /i0  and  /q  are  analytic  in  Qt  except  at 
j  =  --l  provided  fty  is  defined  so  as  to  be  continuous  at  2  —  1.  Without  proper 
choice  of  the  constant  c,  //„  and  /q  will  fail  to  be  analytic  at  —  1.  In  order  to 
render  them  analytic  there,  we  choose  c  so  that 


f  *7i  (■*) 
J  p(s) 


J  9{s) 

-l 


(  Al 

9  V) 


Then  #(  —  1)  —  0,  and  #  may  be  written  in  the  form 


&[s)  =  f  q'{s\~A  ds. 

'  J  2<f(s) 


Provided  iq  is  now  defined  so  as  to  be  continuous  at  —  1,  this  tact  and  computa¬ 
tion  show  that  fi0  and  /q  are  analytic  throughout 


4 


r.*  »<«<%.*  1 1 
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\ 

The  choice  of  uu  and  //,  determine  I)0,  and  from  (4-5)  if  follows  that 
1471  /*•-/»?  ft  al¬ 

l'll  e  relations  (47)  and  (4.1)  an<l  the  definilion  (1.2)  of  R  imply  that 


(4.8) 


h>o(£,A)  =  1  j  A  -  rfhi^rr&y* 


Thus  if  I?,  is  hounded,  it  is  clear  from  (4.8)  that  D„  is  bounded  atvay  from  zero 
for  |  A  j>  .V  and  for  r  £  .  If  is  unbounded,  we  assume 

(iii)  D0  is  bounded  away  from  zero  for  £&t  and  for  j  Aj  .>  N. 

We  observe  directly  from  (4.1)  that  //  -D<(.  Therefore,  we  may  make  the 
change  of  variable 

(4.9)  s -I>0^ 


to  remove  the  first  derivative  term  in  equation  (4.2).  It  then  takes  the  form 


(4.10  a) 

£)-[*’(--■- 

1)  +  A</,(r)  +  T(z,  A) 

]f  =  0, 

where 

(4.10  b) 

+ if. 

+  T-J  4.3/W) 
D0  1  4  Uv 

2  Ido 
2l\ 

The  function  T  is  analytic  in  Si  as  R,  /<0,  and  /i,  are  all  analytic  there.  Further, 
T  is  analytic  in  A  for  J  A  |  >  .V ;  hence,  we  rnay  write 

(4- 1  f )  r(c,A)-V/;(2)A"f. 

o 

We  call  the  differential  equation  (4.10a)  "the  second  approximating  equation" 
for  the  given  equation  (2.5).  It  is  important  to  note  the  way  the  constants  r;, 
which  were  introduced  in  (J.l),  enter  into  the  functions  The  cf  appear  in  T 
wherever  R  and  R'  do.  From  (4-4)  and  (4-3 a),  we  observe  t?iat  aside  fro?rt  the 
leading  term  R  in  the  right  member  of  (4.10b),  R  and  P!  have  coefficients  of 
order  A  1  or  smaller  in  this  right  member.  It  follows  that  the  constant  ct  occurs 
linearly  in  always  with  a  coefficient  of  +1,  and  is  absent  from  every  /,  with 
;<»■ 

5.  The  related  equation.  We  are  now  in  a  position  to  construct  a  dif¬ 
ferential  equation  whose  solutions  art:  known  and  whose  coefficients  are  identical 
with  those  of  the  given  equation  up  to  terms  of  order  A  K  Inasmuch  as  this 
construction  is  already  well  known  [5,  (?],  we  present  only  the  results. 

With  C  standing  for  any  solution  of  equation  (4.10a),  we  set 


(5-la) 

where 

(5.1b) 


r-Dr»(^C+^)f 


A(z,  A)  =  £  «,(*)/-'. 


B'(-'7-)-l6y(:)A--, 

0 
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and 

(5.1  c) 


D,  (.-.A) 


A  A'  -j 


n  s 
;.s 


«  ,j  .  IV 

As  ‘  ’  /.* 


and  choose  the  functions  af  and  ft,  so  that  a-  satisfies  a  differential  equation  of 
the  form 

(5.2)  3S&|«-o. 

where  Q{z,  X)  is  bounded  for  z^fy  and  for  |  A  ]  >  .V.  The  proper  choices  of  the 
a/ s  and  ft/ s  arc: 

<*oW  *  1  •  M*)  *=  4)  j 

1 

r 

(Jt(r)  -=  0,  ftt (z)  —  —  j  [(fo  —  f,)  —  2ql  b0  —  17,  ft0]  — jry  , 


(5.3)  and  for/ —  2, 


^  /"  ('//-*  ~  v-*  ^4. 

1 

.1 

M-1  «  ’  J"  jt  (ft-A+a  -  (,•-*)  «*  -  (2?i  ft,'-i  +  <?i  ft/- 1)  - 

S  -r  f/-* -  s  ft*)  '  j  ' 

In  order  to  render  the  «/s  and  ft/s  analytic  in  £%,  the  constants  r,  occurring  in 
their  defining  expressions  must  be  properly  chosen.  Tliis  may  be  done  recursively 
since  c,-  occurs  linearly  in  f*  and  is  absent  from  each  with  /<«’.  The  correct 
determinations  arc 

1 

f  -  2,h  K  ~  ll  ft0]  J  *  , 


L 

*0 ~y  f  (?* ~ ^0 ) 

ff’.V  ,*  r 

2*'  c'~Y_l  l'h 

(5.4)  and  for  /si  2 

1  -  1 

c,=*yf  |2  (?/-*+* 

-  *t-k)  *»“  (2?I  ft,'-! 

-sW* 

2<p 


where 


(,  =  tj  -  -  Cy  and  y 


5  .  -1 

(/£)  ■ 


We  now  observe  that  the  division  by  D»  in  (S.i  a)  is  legitimate  if  is  bounded,, 
since  l\(z,  co)  —  i.  If  &  is  unbounded,  we  assume 


I 

I 

f 
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(i v)  1),  (x,  X)  is  bounded  ana y  from  zero  for  :  (z  3,  and  (or  j  A|  >  A’. 

The  function  [)(:,  ?.)  which  appears  in  equation  (5.2)  is  a  function  whose  specific 
form  may  he  computet!,  although  it  is  of  no  interest  to  us.  It  involves  the  several 
functions  A  B,  Q,  R,qn,  and  r;,:  It  is,  however,  important  to  note  that  i)(z.  A) 
is  analytic  in  7.  (or  '  7  j  >  A*  and  is  analytic  in  z  lor  zf  3,.  We  henceforward 
refer  to  the  differential  equation  (5-2)  as  the  related  equation.  In  Part  11  we 
shall  single  out  certain  of  its  solutions  and  describe  their  behavior  for  2> 
and  |7|  >  Ar. 

Part  II:  Solutions  of  the  Related  Equation 

6.  Introduction,  The  construction  just  green  of  the  related  equation  enables 
one  to  know  the  behavior  of  its  solutions  only  as  well  as  one  knows  the  behavior 
of  the  Weber  functions  wliich  are  solutions  of  the  first  approximating  equation 
(3.2).  These  are  Weber  functions  of  large  complex  order  and  unrestricted  argu¬ 
ment,  e.g.,  flj(l_n(|'2/! :).  Their  behavior  is  not  yet  completely  known.  Erdei.yj, 
Kennedy  &  McGregor  have  derived  the  results  most  useful  to  us  [.31.  They 
have  established  approximations  to  the  various  functions 

(6.1)  y„(x.r)  ="  v  =  (-  1)", 

with  an  error  involving  v  *.  We  strengthen  their  results  by  showing  that  their 
asymptotic  representations  hold  uniformly  in  aig  v  over  sectors  of  the  r-plane. 
We  also  extend  their  representations  to  asymptotic  expansions,  the  terms  of 
which  are  determined  recursively  by  a  quadrature  at  each  step. 

We  obtain  the  behavior  of  solutions  of  equation  (}.2)  by  letting  *  *=  j  ^  z 

in  the  expansions  for  the  functions  (6.1).  Wc  then  determine  the  asymptotic 
expansions  over  3.  of  the  solutions  of  the  related  equation.  These  are  explicit 
up  to  terms  invoking  v~l.  While  out  results  arc  limited  by  the  lack  of  precise 
information  on  the  structure  of  Weber  functions  of  large  complex  order,  ap¬ 
proximations  to  terms  involving  r~l  are  usually  adequate  in  applications.  When 
they  are  not,  numerical  methods  for  the  calculation  of  the  terms  of  the  asymptotic 
expansions  given  in  §  8  may  be  employed  to  give  more  precise  approximations. 

7.  Domains  in  the  a: -plane.  The  determination  of  the-  regions  of  validity  for 
the  approximations  to  solutions  of  the  given  equation,  which  are  derived  in 
Part  III,  depends  upon  the  character  of  a  function  0(x).  It  is  defined  by  analytic 
continuation  from  its  positive  values  on  (1,  oo),  which  are  given  by 

(7. 1)  0(x)  ~  /  {l2  -  i)4  dt  =  /  <p(t)  dt . 

i  x 

It  is  shown  in  \3,  pp.  469,  470]  that  0  is  a  schlicht  mapping  of  a  Ricmann  sur¬ 
face  X  over  the  x-plane  with  branch  points  at  T  \  onto  a  Riemann  surface 
over  the  0-plane  whose  branch  points  are  the  images  of  the  branch  points  in  X, 
Clearly, 

(7.2)  0(.t)  -  i  {*(*»-  I)4  -  In  [x  -f-  (t*~1)»]}: 
and  for  large  |xj, 

(7.3)  &{x)  «  ±  .Vs  [ 5  +€>{*-*)]■ 

iou 
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The  WdxT  functions  (6.1)  an*  entire  function-  of  their  aigurnent  and  their 
order.  Hence,  in  discussing  them  we  may  assume  that 


(7.4)  Ttie  .r-planr  is  cut  from  -  1  to  -*-  1  and  J  arg  r|  rr  r,  e>U 


These  assumptions  enable  us  to  avoid  considering  more  than  a  finite  number  oi 
sheets  of  the  logarithmic  Riemann  surface 

Let  X,  denote  the  .v-planc  cut  in  accordance  with  (7.4)  and  the  condition 
jarg(x  !)]  tF.*r.  Tim  corresponding  portion  *},  of  3  consists  of  two  sheets  joined 
„  x'-horgx— Jf  alon£  lhe  ra>'  arg <P  —  Q  (see  Fig.  1).  We  now 

fr - — ~ — * - — —  en)t--zi  (jcfine  certain  regions  of  X  on  which  the  be¬ 

havior  of  the  Weber  functions  (6.1)  is  con¬ 
veniently  describable.  Let  the  index  tit,  which 
appears  in  (6.1),  have  range  0,  if,  ±2,  i3, 
pL - - *-tng4-o  and  let 


-crg$-0 


■w-ctg  f?"i> 


£(*,»')  =  v0(.r). 


x<-r,<tryx-jr 

l!lg.  1.  Th»  »irfa<<!  3i 


-erg 


Let  yM  l>e  the  point  at  infinity  on  3  *n  the 
direction  arg  (f>  —  —•  in :t  — •  arg v,  and  let  f  be 
any  positive  number  small  enough  to  make 
the  following  definition  meaningful.  We  define 
61*.  to  be  the  largest  dosed,  simply  connected 
region  Co  of  X  such  that 


(7.6a)  —  §  7i  -}-  e  %  arg  £  4-  in  n  fs  |  ;r  —  c ,  x  £  61 , 

(7.6b)  |x  +  l[  r ,  *6®. 


(7.6c)  |0(*)|£«,  x  £ 6),  except  that  x  —  1  and  a  neighborhood  of  x  —  S  may 
be  in  (#. 


(7.6d)  Each  point  x  in  61  may  be  joint'd  to  the  image  of  yK  by  an  analytic, 
simple  curve  lying  in  C4  and  on  which  5 H£  is  monotone  (monotone  in¬ 
creasing  from  the  image  of  j>„  to  .v  if  m  is  odd,  monotone  decreasing  if  m 
is  even). 


The  regions  61..,  are  illustrated  in  [3.  p.  482  .  In  general,  except  for  sectors  with 
vertices  at  a  —  ±  I  and  widths  2c  in  arg(.v±l)  and  the  interior  of  the  circle 
at  x  =-  —  1  deleted  by  hypothesis  (7  6b),  each  region  61w  covers  the  x-plane. 
Exception  occurs,  for  example,  when  in  0  and  arg  >■  =-(.). 

We  denote  by  65*  tlie  reflection  of  61m  in  the  origin;  that  is, 


(7-7)  61*  -{.vjc  *' if  arg.v^O,  r* ‘ x £ 61  w  if  arg.v  s:  o}. 

Now 

(7.8)  0(e‘-*'x)  =0(x)  f 


so  that  is  precisely  ©„  plus  all  but  a  small  sector  with  vertex  at 

x  -  —  1  and  width  2r  in  arg(.v—  1)  of  the  neighborhood  of  x  —  -  f  not  in  0fw. 
The  definitions  of  61w  and  61®  are  essentially  the  same  as  those  given  in  [3 
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8.  Asymptotic  expansions  of  Weber  functions.  In  this  section  wo  dcscriU* 
the  asymptotic  U-havior  of  the  Welnr  functions  (6.1)  over  the  regions  and 
•"*.  I'li'  se  functions  ire  solutions  of  the  transformed  Weber  equation 


(«■! 


2  -(|S-i)v^o. 


Since  they  are  entire  functions  of  v,  vM  —  yt  if  rn  I  (mod  -4).  fhe  Wronskians 
of  certain  pairs  of  these  solutions  ate  /,  p.  42  4 

-d.-reei""’' 

n,  — o  »--f  t  ,  ‘ 


(8.2a) 

*(vm.  >•«**.  a)  •= 

(8.2b) 

’  y«ii ,  •*) 

where 

(8.2c) 

(!•  *)  •-■= 

and 

(8. 2d) 

a  == 

<ip  <1/ 


We  henceforth  adopt  the  convention  that  in  formulas  in  which  the  double  signs 
i  or  T  appear,  all  upper  signs  or  all  lower  signs  arc  always  to  be  used.  It  follows 
from  (8-2)  that  any  two  of  the  four  distinct  solutions  ym  are  linearly  independent 
except  for  v  —  0  and  certain  other  integral  v  alues  of  v.  If  »•  is  not  zero,  and  we 
always  consider  jvj  to  be  large,  the  solutions  yn  and  ym±l  are  linearly  independent 
without  exception. 

The  asymptotic  representations  for  the  functions  ym (.v,  v)  involve  the  Airy 
function  Ai (/),  which  is  defined  by  the  relation 


in  which  A’j  is  a  modified  Bessel  function  of  the  th.rd  kind  [3,  p.  46)].  For  each 
admissible  ft the  function 


(8.4) 


(0 


C>! 


has  a  simple  asymptotic  expansion  for  ,v  C:  W,., .  Thi  s  is 

(8A.I  A;(r|P)*2-‘3-lr-»fl(r)r4«  +  ° (lf 1 ' 

»«  -0.  if.  ±2,  ±3- 


The  asymptotic  series  for  dXjdx  for  x  £  (V>m  may  be  obtained  from  (8,5'  by 
differentiation.  Ekdei.yi,  Kennedy  &  McGregor  have  shown  that  the  func¬ 
tions  (8-4 )  are  asymptotic  representations  for  the  W.-b<  -  functions  (6.1}  [3,  p.  479j. 
These  representations  may  be  extended  to  asymptot  c  expansions  by  using  the 
algorithm  given  by  Lancer  in  [5]  or  similar  algorithms.  We  use  Lancer's. 
In  the  notation  of  [5j,  it  enables  us  to  construct  certain  functions  and  ( i f 

1  The  form u Lis  in  [/]  apjuar  to  contain  misprints. 
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which  arc  involved  in  the  coefficients  of  the  asymptotic  expansions  for  the  Weber 
functions  (6.1).  With  reference  to  equation  (8.1)  tin-  defining  expressions  of  a, 
and  /?;  specialize.  7  hey  Ix-come: 


(H.6a) 

and  for  /2>  2, 


*«(*)  ~  1 , 
a,(x)  ~  0, 


i 


fiiix)  ~  d* 


d  I , 


V  f  W'-tV) +k(i)fiHt(t))dt. 

(8.6  b)  *  „ 

ftW  -  ^  /{«"«  +A0M)  +2^-.«3  -*W/-.W}  41/)  - 

I 

wherein 


(8.6c)  k(x)  =  and  ^(.r)  =  01  (*)?-*(*). 

Since  cq— A  — 0,  a5/+1  — — 0  for  all  /SrO.  The  following  order  relations  are 
easily  established,  They  are  important  in  the  extension  of  the  asymptotic 
representations  in  [.3]  to  asymptotic  expansions.  For  |  a  |  >  A7  and  ;2>  0, 


(8.7)  «*/W“0(l).  and  &,(*) -0  (*">). 


The  asymptotic  expansions  given  below  extend  the  results  in  [3j.  We  omit 
their  derivation  to  avoid  repetition  of  work  in  [3].  The  a/s  and  /3/s  are  suf¬ 
ficiently  small  at  x  —  oc  to  guarantee  uniform  convergence  of  the  integrals  which 
enter  into  the  derivation.  This  is  the  point  in  the  extension  of  the  analysis  in 
[3]  that  requires  careful  attention.  The  bounds  for  the  error  terms  in  our 
expansions  are  uniform  with  respect  to  arg  v,  if  arg  v  lies  on  a  bounded  interval. 
This  also  applies  to  the  results  in  [3],  The  reason  for  this  is  that  the  limit  am 
of  integration  in  [3,  equation  (4.8)]  may  be  fixed  provided  arg  v  has  a  range  of 
length  less  than  n.  For  example,  if  |argr|  \  jc  —  e,  aa  may  be  chosen  as 
x—  r0'*  oo;  and  if  jargv  —  £r?|js>  J:r  — r,  a0  may  be  chosen  as  x  —  c~^ni-  co. 
For  each  such  range  of  arg  v,  a  bound  on  the  integral  in  [3.  equation  (4.8)]  may 
be  found  which  is  independent  of  argr.  It  follows  that  for  arg  v  bounded,  bounds 
for  the  error  terms  independent  of  arg  v  can  be  found.  This  reasoning  does  not 
imply,  however,  that  in  the  x-plane  the  regions  of  validity  for  the  asymptotic 
expansions  below  are  independent  of  argr.  They  are  not,  although  in  the 
neighborhood  of  infinity  the  boundaries  of  the  regions  of  validity  in  the  ]vx- 
plane  are  asymptotically  independent  of  arg  r. 


(8.8  a) 


y«(*.  »’)  =Cm[/J?(.v,  v)] -i  x 


X  j/1  *,(x, v)  X„ [x, v)  [1  -}-0(»-“/_1)]  -f 

for  x  (E  ©„ ;  and 


(8.8b) 


x[4%v,r)A'n(x.v)4-i?,JV)  +0 (v~'~')}, 


dx 
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HI 


(8  9  a) 
(8.9  b) 

and 

(8.9c) 


0.VO 


Cm  us  »d *  f"  I  «  »>i‘ 


.■i*(r,  r)  v  1  ll*(x,  y) 

Am'{. r,r)-(**-i)  B*[x,r)  k[x)  ft*(*  v)r  !  A*(x, *■}  +  *'*  //*'(*, r)j‘ 

A  *(.v,  »•)  » 'v  a,  (a)  v-.  /*•(*,  r) »£/?<(*)  v~h 

C-  0 

The  functions  a,.,/),,  and  A  arc  defined  by  the  relations  (8.6),  the  functions 
A'„(.e,  v)  are  descrilx'd  by  the  formulas  (8.4)  and  (8.5),  and  a~  {—  1)”. 

One  uses  the  identity 

(8-50)  ym  (x,v)  =ym±t{^!'ix,y) 

to  derive  the  form  of  ym(x,  v)  when  a  C  03®^,.  Taken  together  the  regions  f*}OT 
and  {»i  =  0,  ±  1,  ±2,  ± })  cover  the  x- plane,  each  point  being  in  two  regions 
having  indices  which  differ  by  1  or  ).  Therefore,  the  expansions  (8.8)  yield  the 
asymptotic  behavior  of  two  linearly  independent  solutions  of  the  differential 
equation  (8.1)  for  each  ,v.  We  observe  that  zeros  of  the  Airy  function  are  not 
excluded  from  the  region  where  the  relation  (8.8b)  is  fulfilled.  This  is  the  reason 
for  the  altered  character  of  the  error  term  in  (8.8b)  as  compared  with  the  error 
term  in  [3,  equation  (9-6)].  for  example.  The  behavior  of  dvjdx  may  be  found 
via  the  relation 

(8-11)  Ld^.^vD,._,(x)-'xD'[x) 

from  the  expansions  (8  8). 

It  is  interesting  to  note  that  if  v  is  real,  all  solutions  of  equation  (8-1)  are 
oscillatory  on  the  interval  —  i  <  x <  1 .  This  occurs  since,  for  the  configuration 
of  x  ana  v  cited,  iff  (£(x,  v))  —0. 

9.  Solutions  of  the  related  equation.  The  first  approximating  equation  (3,2) 
has  solutions 

(9.1)  (m  -0,  ±i,  ±2,  ±3), 

where  D„(z)  is  the  standard  parabolic  cylinder  function  of  order  r  [i],  and  where 

(9-2)  v  —  A  —  c  —  V  A  ~,_1. 

o 

Of  course,  vm  —  v,  if  /s=»i  (mod  4).  The  formula 

(9.3)  “y«(  |  **.*') 

relates  the  solutions  of  equations  ().2)  and  (8.1).  In  view'  of  this  connection, 
the  asymptotic  behavior  of  the  solutions  rm(r, /.)  may  now  be  easily  described 
with  reference  to  the  .r- plane.  Since 

(9.4)  Xss\rZt 

the  transition  from  regions  on  the  x -plane  to  their  images  on  the  r-plane  is 
elementary. 


r'<,W5££i?3SSS?*| 


M2 


Nittinj-vS  h 


i>y  directly  the  theory  in  •>  to  (1.2i,  we  cuulii  have 

described  tin*  Iwtdinvior  of  certain  solution-  of  (  l  2),  other  than  flic  solution*  (0.1), 
on  regions  in  tin*  --plant*  identical  with  the  r«  gion*  and  Of*.  Hmvevei,  the 
slight  simplification  in  tin*  desciiptieu  of  the  regions  of  vaiiditvof  the  asymptotic 
exjwnsions,  which  is  gamed  l>\  this  approach,  i-  offset  f>v  an  increase  m  dif¬ 
ficulty  in  deriving  tin  sc  expansions,  since  equation  (V2<  rs  superficially  more 
complicated  than  equation  (8.1).  For  this  reason  and  to  Sake  fullest  arlvantage 
of  the  work  in  [.?},  we  proceeded  as  we  have  done. 

Wc  compute  the  Wnmskiaiss  of  the  solutions  7.)  using  tin  formulas  (8  21, 
(9.1),  and  (9.4).  For  future  reference,  the  computed  values  are: 


2l.iAc*"a' 


r 


n  r-t-  I 


(9.S) 


To  each  solution  n,B (;,/,)  of  equation  (}.2)  there  corresjronds,  through  the 
relations  (4.1),  (4.0),  and  (5.1),  a  solution  um{z,  ).)  of  the  related  equation.  The 
solutions  re*, (r,  /)  may  be  written  in  terms  of  the  functions  vm(z,  /,)  as  follows: 

A)  ~  £*(;.  /)  vm{:, ).)  ■+•  h't  {s.  /.)  i£{s,  /•)//, 


£0(-*, A)  »(»,«,)-* 


1  /it 


t<  j  !‘,  K  !'«  _  /‘r.  I >o 

}.  >  7*  ‘r  A.  2  A./V) 


s,  ;.)  _  (0.  /»,) - 1 1  ,i ,,,  +  '!  •*  "•  -  |j . 


(<).<>) 

where 

(9-7) 

Since 
(9-8) 

the  linear  independence  of  solutions  tr,  and  try,  depends  upon  that  of  their  cor¬ 
respondents  ty  and  vm, 

For  the  proofs  to  follow  it  is  necessary  to  have  at  hand  the  asymptotic 
behavior,  in  first  approximation,  as  J7.|  -*•  co,  of  the  solutions  u'm[z,  ?.)  of  the 
related  equation  together  with  the  behavior  of  their  first  derivatives.  One  finds 
this  behavior  by  considering  the  relations  (9.6)  in  conjunction  with  the  asymptotic 
expansions  (8.8)  and  (8.5)  and  the  relation  (8.4).  When  j  ,f(.v,  r)|;£  M  and \x—  l  j  < n 

(9-9a)  u,„(z,Af  ■=  Cm 0(1),  and  iv'm (:,/)  =  C„0 {#) . 

When  j £(—  x,  »>) j  and  j.v  -}-  lj  < c. 

(9.9b)  ^'(U)=C.±tO(1).  and  u,-'m(z,  A)  -  Cm±a0(Ai). 

If  i  is  bounded,  |£  (.v,  r)|  >M,  and  |£(-x,j>)|>Af,  then 


(9- 10  a) 
(9.10  b) 
where 
(9-IOc) 


»»M)  =  7)  (A (A- 1), 

hi(:,  7.1  u-,„{z ,  7.)  O(A) , 

A„(z,  7.1  -  id"rc  l •’•■ii «( 


mKmammmmrn 
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If  is  not  boundid,  then  the  functions  and  /;,  defined  by  (9-7}  may  lx*  un 
Itomvied.  In  tins  c.v>e,  tin-  relations  (9-loaj  and  (y.lOb)  aic  invalid.  We  therefore 
introduce  i  Junction  !>(:)  which  is  continuous,  which  is  nonvanishing  for  ,;!>.!/ 
and  equal  to  1  for  .11,  .trid  which  is  the  "largest'1  function  such  that  for  each  m 


'<>.!!  M.1J(;  i)j-  .v;r;  1  ami  >.)  : 


u  , h  j 


•  .v  »  ;-i. 


whin  t  blw  t?C<«  and  when  £  (a.  s*l(  >  .1  /  and  jc{-  a,r)[>jl/.  Thus  when  the 
rvlnti  >ns  (y  1 1)  apply, 

...  /.  u  _ 

and 


(9  12. i) 

and 

(9.121)) 


"  b(z) 


Part  III:  Asymptotic  Expansions  for  Solutions  of  the  Given  Equation 
10.  Final  hypotheses  on  Our  objective  is  to  determine  the  asymptotic 
behavior  throughout  S\  of  a  pair  of  linearly  independent  solutions  of  the  given 
equation  (2.s).  This  will  be  done  in  the  sections  to  follow*.  We  first  link  the 
behavior  of  solutions  of  the  given  equation  to  that  of  solutions  of  the  related 
equation  by  using  the  familiar  method  of  variation  of  parameters,  its  application 
to  the  given  equation  yields  the  integral  equation 


10.1) 


u  ~~  w  -j*  j 


'  '*  j  h-.;  a '/,  {') 

;,W) 


which  is  equivalent  to  the  given  equation  in  the  sense  that  an  analytic  solution 
of  either  one  is  a  solution  of  the  other.  In  equation  (tO.f),  wa  and  ws  may  be 
any  pair  of  linearly  independent  solutions  of  the  related  equation,  w  may  be  any 
solution  of  the  related  equation,  and  may  be  any  point  in  The  kerne! 
in  equation  (10.1)  is,  of  course,  independent  of  the  choice  of  <e0  and  xvb  so  that 
a  unique  solution  u  of  the  given  equation  is  determined  by  specifying  w  and  -:s, 
in  equation  (10.1). 

The  variables  z  and  x  and  the  parameters  ?.  and  r  are  always  considered  to 
fulfill  the  relations  (9.4)  and  (9-2).  Thus,  given  z  and  v  and  x  are  determined. 
In  order  to  keep  the  investigation  to  a  reasonable  length  and  to  avoid  further 
notation?.!  complications,  we  make  the  assumption  that 

(v)  jargAj^.-r. 

The  second  of  the  conditions  (7.4)  is  thereby  fulfilled  when  ]7.i  is  sufficiently 
large.  There  is  no  loss  of  generality  in  assuming  that  the  first  of  these  conditions 
is  also  satisfied.  In  what  follows  we  fix  7.,  once  for  all,  with  j  A|  sufficiently  large 
to  fulfill  all  requirements  placed  upon  it. 

The  following  Jemma  will  be  used  in  the  discussion  of  solutions  of  (10.1). 

Lemma.  // 

a)  -ssf  is  a  dosed,  simply  connected  region  (zahose.  boundary  may  depend  upon 
a  complex  parameter  }.), 


ex*  *  r 


if 


- - 


U-J 


Nicholas  1>,  Ka.*amn<>h 


h)  /(:.  k)  i*  cor.Ur.uoHi  and  tnwnJcd  j or  :  .»/  and  for  \/.\  >.Y, 

C)  /i’(.\  7)  ii  {(.w/ in homs  fur  :  Ccsf,  far  I ,  ,«/  wm!  Aw  A  >.V, 

(1)  /Arrr  i*s  is/  <i  constant  M,  a  point  s  t/,  t*m/  a  set  1 1  red  triable  analytic  curves 
flirting  the  points  of  S»  t<<  \  such  tha‘  for  all  :  >■'- .'/ 

/|A'(r.  /.;.)! !«//'  st. 

//it’ii  the  solution  of  the  integral  equation 

gM)  — 7  / (c,  /)  -r  /  x  f  K(:,t./.)g[t./.)dt  (k  > 0) 

i 

/ws  the  form 

?{-’./)=/(*, A)  + so(i), 

where  0(1)  denotes  a  function  of  z  and  ).  which  is  uniformly  hounded  for  :  -ex'  and 
)or\X\>N. 

The  hypotheses  of  the  Lemma  guarantee  the  uniform  convergence  of  the  series 
which  is  obtained  by  successive  iteration  of  the  integral  equation,  and  from  this 
the  proof  of  the  Lemma  follows. 

For  the  Lemma  to  apply  to  the  integral  equation  (10.1),  it  is  necessary  that 
the  kernel  in  (10.1)  lie  bounded.  We  now  make  assumptions  on  £&t  sufficient  to 
establish  the  boundedness  of  the  kernel.  For  each  /,  let 

(10.2)  9tm  —  {;( x  (z)  £  im  —  0,  ;r  1.  ±2,  ±  3)  • 

In  ttiis  definition  and  in  the  sequel,  when  a  region  with  subscript  larger  than  1 
or  less  than  -  3  appears,  the  region  is  understood  to  be  ti  c  null  set.  Because 
of  the  overlapping  of  the  regions  0)M  and  also  of  the  regions  ($*,  one  can  show 
that  j v/ j. j  (m  “0,  ±1,  ±2).  It  follows  from  the  hypothesis  (v) 

and  the  conditions  (7.41  that  we  lose  no  generality  by  assuming  that  df.yddt. ~ 
and  i?3C.^s.  The  regions  cover  &\  and  each  point  of  %  is  in  at  least  two 
such  regions.  Our  final  assumptions  on  S>:  are: 

(vi)  Corresponding  to  each  region  Mm,  there  is  a  point  zm  such  that  each  point 
~G'~m  way  be  joined  to  zm  by  an  analytic  curve  F  lying  in  and  such  that  on  F, 
iK  £  is  monotone  (monotone  increasing  from  ;m  to  z  if  m  is  odd,  monotone  decreasing 
from  sm  to  x  if  m  is  even). 

(vii)  The  integrals 

J  ~  b'V) 

r 

are  uniformly  bounded  for  ::  £  and  for  \  A  j  >  A7. 

The  function  .0(z,  A)  first  appears  in  equation  (5.2).  The  function  b (tj  is  defined 
by  the  conditions  (9.11).  Of  course,  if  £?,  is  bounded,  hypothesis  (vii)  is  auto¬ 
matically  fulfilled. 

There  are  certain  shadow  zones  which  must  be  excluded  from  the  regions 
in  the  discussion  of  subdominant  solutions  of  the  given  equation.  Their  images 
on  the  .r-plane  adjoin  the  segment  --  1<  v<0.  Their  si  re  and  presence  depends 


1 


atititr 
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upon  ,irg  /  They  arc  defined  as  follows. 

-{-!-•  •*,  ■  « ("  £)  >  A'.  \i(  t  {:)}  j  >  A*. 

(10.1)  and  the  paths  joining  :  to  ail  jxtss  through 

the  ncighlHirhiKKl  of  where  \;(  -v)  |  <  A"  j 

The  shadow  /ones  are  excluded  for  the  teason  that  the  subdominant  solution*; 
of  the  given  equation  undergo  a  change  in  their  asymptotic  forms  as  one  crosses 
into  them.  The  presence  of  tiles'.*  shadow  zones  (see  Fig.  2  for  an  example)  is  a 
phenomenon  which  has  no  counterpart  in  the  discussion  of  asymptotic  solutions 
of  differential  equations  having  just  one 
turning  point. 

As  we  have  remarked,  we  wish  to  de¬ 
termine  the  asymptotic  structure  through¬ 
out  of  a  pair  of  linearly  independent 
solutions  of  the  given  equation.  We  first 
prove  that,  for  each  in,  there  exists  a 
solution  um  which  is  asymptotically  re¬ 
presented  by  wm  in  the  region 
We  then  derive  connection  formulas 
among  the  solutions  um  from  which  the  behavior  of  the  solutions  u„  may  he 
derived  for  all  z£@t.  Lastly,  we  discuss  certain  characteristic  solutions,  which 
may  well  Ik*  of  special  interest  in  applications  of  the  theory*. 

11.  Subdominant  solutions  when  v)|  > N.  For  each  region  .#m  — 

there  is  a  solution  of  the  given  equation  which  is  subdorninant,  exponentially 
small,  for  :  f  —  Sfn)  and  Si(c7£)>0.  This  solution  is  unique  up  to  a  constant 
factor.  We  now  establish  its  existence  and  give  its  asymptotic  properties.  We 
assume  for  the  moment  that  z~  {8im  —  S?n)  . A.  l’n  equation  (10.1)  let 

wa—wm  and  wb*swm+i ,  let  tc*  —  um ,  and  let  the  path  of  integration  be  a  curve  V 
which  originates  at  z„.  For  brevity,  we  adopt  the  notation  r  v)  and 

{  —£(x(z),v).  If  the  curve  /"should  intersect  one  or  both  of  the  regions  where 
|r(±x)i  :SA’,  then  $(cr£)i£0.  Thus,  if  we  replace  the  subarcs  of  F  on  which 
j  r  (i  a) j  gj  X  bv  arcs  on  which  j  r  ( +  x)\  —  Ar,  and  if  we  note  that  on  the  remaining 
portion  off,  3i  (a  £)  $  (a  r) ,  then  we  see  that  on  the  (new)  path  of  integration 

€.\p[2o(£  —  t)]  is  bounded. 

The  choices  just  made  determine  a  solution  t.i„  of  equations  (10.1)  and  (2.5). 
If  we  adopt  the  abbreviations 
{  U„(z)  ^A;!b(z)uJz,A), 

(1M  u;(z)-/l  A). 


a>y  9 — tx 


Fltf.  2.  /  ihailuw  tow. 


where  Am[z,  A)  and  h{z)  are  defined  by  (9.10c)  and  (9.13).  we  can  rewrite  the 
integral  equation  (10.1)  in  the  form 


(11.2)  U„  (z) «  Wm  (a)  --/  {  lYm  (z)  ,  (0  -  Wm  {<)  Wm,  x  (a)  a* 


rnO[i)iJ(l)UM(r,  , 


In  writing  (It. 2),  we  have  used  the  evaluations  (9.8)  and  (9-5)  of  ^"(uy, ,  u\n  .  r}. 
The  relations  (9.10),  (9  12)  and  (11.1)  imply  that  for  |  v)j  >  A"  and 


i 


i 
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:  (.V„  '/j  Ix.tli  llw  and  HL. ,  an-  Wi  haw  already  ub'i  rved 

lli.it  t'.vp  2o(;  r!  i-  i<iundid  on  }'  All  function'  iti  equation  (11.2  .tii  tun- 
tumouv  In  urtiK'  of  these  facts  ami  (m)  "f  $  10.  tin*  hy|*>the'»  s  of 

tin*  Lemma  an  fulfilled  with  irsixct  t«»  equation  r  i  1  2).  Tints  for  Ym) 

V„  ]‘(  j  »(:),r)i  \  and  ).  A'. 


We  ran  show  In  analysis  analogous  In  that  above  that  this  relation  a  No  holds 
for  ; .#m  ,  In  trims  of  um  aiul  the  result  (1 1  A)  takes  the  form 


».<U)  -'i$T  (M!>-V)- 


(11.4)  hold'  for  all  ^(Ji',,  such  that 


Since  (.#m  .9’m)  ,  >Jn  , ,  (11.4)  hold'  for  all  ^„)  such  that 

|S(±-v(c),r)j>X 

Wc  can  iK  the  estimates  f  1.4)  and  (0-12)  to  obtain  the  the  form  of  u^(;,  /.) 
directly  from  the  equation  found  by  differentiation  of  (10.1).  It  is 


)>  :  v  ■  •  lM  (;,/.)  O  (l ) 


where  the  prime  denotes  differentiation  with  respect  to  ;. 

12.  The  solutions  «m  when  v)|fS.V.  There  are  two  cases  to  he  con¬ 

sidered:  |£ {x,  >')(  *>  N  and  j£  (  —  x.  v)|  <  A’.  For  the  moment  we  suppose  that 
j  *  j  -_;I  A*.  We  recall  that  the  solution  of  the  given  equation  and  equation 
(10. t)  is  that  solution  determined  by  the  choices  :r— trM  and  in  (10.1), 

the  path  of  integration  being  a  curve  f.  We  divide  f  in  two  parts  /j  and  l',t 
on  which  |  r j  i?: N  and  |t|<A’,  respectively.  As  was  done  in  §11,  we  replace 
any  part  of  1]  on  which  Jr(  ,i)j<AT  by  an  arc  on  which  |  r ( — .v) j  -  A'.  We  now 
write  the  integral  equation  for  u.„  in  the  form 


(12.1)  um(:)  +  /  K(z,  t,  X)  um(t)  it  +  f  K(z,  I,  X)  ujt)  dt, 

!\  r, 

where 

A'/-  f  Jj  _  “'inW  1  (0  ~  N'm(<)  Il'miifl)  Q(t.X) 

“  »>*.« Wi,. n  '  V"  ' 


the  upper  or  lower  signs  being  used  according  as  /J  j  or  TJ  C ,. 

For  the  sake  of  argument  let.  /,  '  On  i] ,  the  behavior  of  uln  is  then 

described  by  formula  (11.4).  Using  this  result  and  the  relations  (Q.9a),  (9.10), 
and  (9.12),  we  find  that 

/*■(.-.  U)  »„</)■«-  Cic,.,  f  0(1)  t 

J  I  /  I  .  I  r  />sp)  Aw'-y 


f[0{i)+C  3"rO( t)j^)  dt. 


Hut  l)i  (or)  A  0  on  !\ ,  except  perhaps  for  an  arc  on  which  flf  (m)  A:  -  ■  A’.  Therefore 
by  hypothesis  (vii), 

(12.2)  j  K(:,  t.  X)  X)  dt  ^  (?n<;(20. 


89SMtsaeiB@i**^ mss  w ,»rv»wrjn«5^.e  r  icv^rw 
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1  mm  tin  I'tiMuic  mt!  the  <ird<  r  relation  I «/#  <>(/.  wc  dedin  >• 


M2.li 


/  h'C.l.  /.), It  /  *  ‘  /  *  '*  ‘  0{ ! ) 


I  he  estimate  (12  2)  it-  to  nuritf  equation  (12.1)  a* 

" "  V"  •  /  '  ‘<){))  .  f  A(;,  /,  /)  "m(M  <ll. 

Hit*  i  flat  ion  (12  I)  shows  that  tin*  hvjxjtheses  of  t  lie  Lemma  arc  fulfilled  uni. 
t rsfK'Ct  to  this  equation;  hence, 

(12.4)  «„(*',/)-  ^(:,/.|  f  (|4(a,»-)!  A'). 

Since  tlu*  order  relati  in  (0.9a)  is  independent  of  arg(r~  )  >*/A),  this  result  is  valid 
in  a  full  neighborhood  of  x  — *  1,  or  of  z  —  J  17  A ,  w  here  jf|£»A'  and  not  just  in  that 
portion  of  Jtm  where  |c|  S  A’. 

We  now*  consider  the  case  |£(  ~.v)j  g  A\  This  situation  is  more  complicated. 
If,  for  such  31(a$)<  —  A\  then  the  estimat  •  (12.2)  no  longer  holds,  and  therefore 
we  cannot  determine  tfie  behavior  of  »<„  when  |£(  — a))5CA'  as  we  did  above  for 
1  >  I  §2  A’.  We  resolve  this  difficulty  by  showing  that  it  is  only  an  apparent  one. 
Let  us  momentarily  suppose  that  m~  T2  and  3i  (<;  £)  -  N.  Then  3i  (a  $)  >  A’ 

in  the  neighborhood  of  Jfm.{  where  j^(  By  the  argument  given  in 

the  next  paragraph  we  are  able  to  determine  the  behavior  of  u„.x  in  a  full 
neighborhood  of :  —  —  |  r/A,  oi  of  v  -  =  —  1,  where  j  £(  r)|  %  A’.  We  are  also  able 

to  determine  the  behavior  of  um  ,  in  such  a  neighborhood.  Similar  considerations 
apply  to  values  of  in  other  than  ±2.  Thus,  we  shall  always  be  able  to  determine 
the  behavior  of  two  linearly  independent  solutions  of  the  given  equation  in  a 
full  neighborhood  of  :  —  —  Jr/A  where  l£(  -v)j'>A’. 

As  promised  above,  we  now  obtain  the  form  of  /r,  if  i£(  v)i *V  and  if, 
for  such  z.  lH(of)S  -  A*.  To  do  this  we  replace  um {:)  and  T.-m n , (r)  in  (12.1)  by 
u’nb2(t ' 1  ‘ :)  and  trm , ,  (c  z),  respectively,  and  similarly  replace  wn(!)  and 
a*m  i.,(/)  for  r  such  that  j  r  (--  *)[  <  .V.  In  virtue  of  the  identity  (<S.  10),  the  analysis 
used  in  the  derivation  of  (12.4)  may  be  easily  modified  to  yield  the  conclusion  that 

(12.5)  «„(U)  -  +  Cy^']  . 

when  |£(  x,  »*)j g  A’  provided  that  for  such  z,  31  (a £) Sr  — ■  A’.  In  (12.5)  the  upper 
signs  are  to  be  used  if  m<0,  lower  signs  if  /«>(),  and  either  if  w  o.  As  is  the 
case  when  j£|<A\  we  need  not  impose  the  restriction  here  that  «  lie  in 
We  need  only  suppose  that  the  portion  of  the  path  of  integration  on  which 
j£(-  ,v)|<;  Ar  is  rectifiable.  Tims  the  result  (12.5)  is  valid  in  a  full  neighborhood 
of  z  |  >7 A  where  J  i  1  —  .v)[  rg  A'. 

We  ran  use  the  estimates  of  um  obtained  thus  far  to  find  the  form  of  A) 
when  |£j<IAr  or  |f(  —  jr)jg.V  directly  from  the  equation  obtained  by  differentia¬ 
tion  of  (!0.1).  The  results  are: 

(116a)  «»M)  -aCM-r  (|iU..')|S.V), 

(1161))  ,(r"*a  *•  (l«  '.Cl-'.V). 

in  neighborhoods  of  z  =-  j  r/A  and  :  —  -}y  A,  respectively. 


Wc  sumnuri/v  utir  results  in  i  theorem: 


Theorem  1.  I'nder  the  hypotheses  (i)  through  (vii)  and  corresponding  to  ear/; 
region  ,  the  givrtt  equation  (2. 5)  has  a  solution  nn  with  the  following  asymptotic 
behavior  when  |  A  |  >  .V. 

a)  li'A/n  i t  (•*„ •—  y„)  mid  je(d::-v,  r)J  >  .V,  um  n  exponential!',  large  or  small 
according  as  3{(»7i)*  0  rr  4{ {o z) > O  and  is  described  bv  the  formulas  (11.4)  and 

(M.5). 

1>)  H7;e«  |;  -  |  vjX j  *£  /  and  t(x,  »)|  g X  or  j:  -r  1 1//.|  <C  r.  and  jc(  — a.  >•)(  gj X, 
um  is  oscillators,'  in  character  and  is  described  by  the  relations  (12.4),  (12.6a)  and 
(12.5).  (12.6b),  respectively,  provided  (»•//.  and  —  Jr//,  lie  in  ,it„  and  'M (//£)  el  —  A’. 

c)  If  v  is  real,  then  um  is  bounded  and  oscillatory  for  ~  1  ^  .r  ^  1 ,  n;  —  0,  £- 1 ,  i  2. 

The  regions  3fm  and  S/m ,  the  variables  x  and  £,  and  the  parameter  v  and  a 
arc  defined  by  (10.2),  (10.1)  (9-4),  {?■>),  (9-2),  and  (6.1),  respectively.  The 
functions  and  constants  appearing  in  the  relations  (11.4),  (11.5),  (12.4),  (12.5) 
and.  (12.6)  are  defined  in  Parts  I  and  II.  'flic  hounds  for  the  error  terms  in  all 
of  the  estimates  for  solutions  of  the  given  equation  which  have  been  derived 
arc  uniform  in  arg  a  for  jarg/.jSsrr.  This  is  so  for  the  same  reasons  as  those 
cited  in  §  £  for  the  uniform  boundedness  of  the  error  terms  in  the  expansions 
(8.8a)  and  (8.8b)  with  respect  to  arg  >•. 

J3.  Connection  formulas.  It  remains  to  give  a  description  of  for  solutions 
/„(■?,/)  when  | s? ( ±  .v) j  >  A'  and  whc.i  z  is  not  confined  to  the  ipediic  region 

as  assumed  in  Theorem  1.  The  dependence  relation  among  three  solutions  :ih 
uM,  and  «w+,  of  the  given  equation  may  be  written  in  the  form 


(13.1)  Ml  A,  „um  +  B,  „um+1. 

The  coefficients  are,  of  course,  given  by  the  formulas 

(u  H  A  _  _f(«|.»0  +  l,£)  n  _  tlt, 2) 

K  ~  ir\um. »m+!,r)  • 

Since  the  given  equation  lacks  a  first  derivative  term,  the  Wronskians  in 
( 1 3-2)  are  all  constants  and  may  be  evaluated  at  z  —  ]  rj?. ,  where  x  --  !  and  £  =  0. 
If  each  of  the  solutions  u,,  um,  arid  umi  j  fulfills  the  relation  (12.4),  then  we  may 
conclude  that 

*>*.  »/)  =■-  ^(wA.av)  +  C*C,0(^-— *). 


where  k,j=l,m,  or  Thus  by  relation  (9.8), 


(13-3) 


Al,tn  —  d'-rn+l  U\''  )> 

V  V  m.  ‘  ,H  i) 

Ji  —  h'm-  ,  ('  r~  A  a  n  1 1  ()  I  J  -  n-S  S\ 


Since  vk  —  Vj  if  As=;  (mod  4),  the  relations  (Q.5)  always  suffice  to  give  the  values 
of  the  Wronskians  in  (13-3)-  The  constants  Cy  are  defined  by  (8.9a). 

The  indices  l  and  m  for  which  the  formulas  (13.1)  and  ( 1 3 . 3 )  .hold  depend 
upon  argr;  hence,  the  pairs  of  linearly  independent  solutions  of  the  given  equa¬ 
tion  whose  behavior  can  be  determined  over  all  of  depend  upon  arg  »•.  The 


j.  Him  rr«irrm«wiir  rr  -i - 1 — rm — rmn — 
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selection  of  these  fundamental  juit*  of  solutions  may  lx.-  made  using  Theorem  I 
arid  the  connection  formulas  (1J.U.  Wc  omit  the  detail*  of  the  selection  and 
-t.iti:  ord\  the  trstllts 

Theorem  2.  Tor  each  r  line  of  /  suck  that  jarg  ?.\  1$  .1,  the  given  equation  (2.S) 
has  a  pa\r  of  linearly  independent  solutions  ichose  asymptotic  behavior  over  all  cf 
'Jj  may  be  found  from  Theorem  1  and  the  connection  formulas  (It  t);  namely. 


0  arg  r  <  ,n  | 

«  «,»i 

rr  S.  argr^  n  ~r  i  j 

Ike  fundamental  pair  is  • 

H-t.  “o 

rr  <  arg  r  0 

II  , ,  Hj 

—  .7  —  f.  arg  v  g  —  n 

M0.W,. 

14.  Characteristic  solution.1;.  If  J{irr~  l)=A'.  a  non-negative  integer, 
>’«(*.»’)  =  i  //*  (}  1 A  -f  2x), 

where  y„(r,  r)  is  a  solution  (6.1)  of  the  transformed  Weber  equa'ici  (8.1)  and 
//,  (a)  is  the  Hennite  polynomial  of  order  k.  The  Hermite  polynomials  are  even 
or  odd  according  to  the  parity  of  k  so  that  when  l  (or-  1)  k, 

----  (-  i)*yn±2(.v,r). 

These  facts  arc  reflected  in  the  behavior  of  the  solutions  of  the  given  equation. 
If  j arg / j < f , 

^(Ko.m**)  =  »:'*f-»rc  ,:t,[0(;.  M-‘)  0  (/"-’)]. 

where  the  upper  or  lower  signs  are  to  be  us.d  according  as  arg  r<  0  or  arg  »•>(*. 
When  r~2k  -f-1,  T/'(v0,  y_ ,)  —  0.  Thus,  as  shown  by  Lanouk  [4],  for  certain 
values  of  /  the  Wronskian  H'(u0,  u  _.)  vanishes.  These  characteristic  values  are 
countable  in  number  and  are  given  by  (he  formula 

2k  1  r  OUr  ')• 

At  the  ciiaracteristic  values, 

«0M*)  ==(-i)*«±*(U*)  [i  +o(;.4- -•)]. 

where  the  i  oeror  lower  signs  are  to  be  used  according  as  arg  (>'(/.*))  is  negative 
or  positiw.  ih  <  i  n,  '  I  a  ''haracteristic  solution  of  the  given  equation. 
We  can  observe  u.e  u  i.  -  ot  a  characteristic  solution  u0(z, for 

real  7.  by  referring  to  Theorem  1 . 

Since  |$f(A*$)j  <»V  for  values  of  r  such  that  —  ts£  5*  1  an  1 1 A  - 1  <  .V  |  ?.k\  *, 
a  characteristic  solution  is  oscillatory  for  :  in  a  strip  of  width  Op  )  centered 
on  the  interval  [—  \,  fi.  in  certain  regions  outside  of  such  a  strip,  a  characteristic 
solution  is  subdominant  and  decreases  exponentially  to  zero  as  |r|->-<x>.  These 
regions  are  sectors  of  width  \ .7 -  c  with  vertices  at  z  —  ±\/.Jv(2.t\  and  which 
arc  centered  on  the  real  axis.  In  particular,  a  characteristic  solution  is  uniformly 
bounded  for  real  z.  Only  characteristic  solutions  and  their  multiples  have  tins 
property. 
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